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ABSTRACT 


The  electnm  propagator  theq^  is  presented  with  scmiewhat  of  a  historical  perspective  and 
the  woridng  equations  are  developed  with  the  aim  to  take  advantage  of  molecular  pcwt  group 
symmetty.  A  new  electnm  propagator  code,  the  vectorized  electnm  propagator  program  (VEP), 
is  intio&ioed  without  full  details  about  its  structure  and  cqmbilities  (such  details  are  bring 
published  elsewhere).  Applications  to  dm  (UV)  photoelectron  q)ectra  of  some  donor-accq>tor 
complexes  of  borane  with  carbon  monoxide  and  water  are  presented  at  the  level  of  second-order 
dieoiy  as  an  illustration  of  the  theoiy  and  die  VEP  code. 


94-15930  J, 

4$*'  94  5  26  1  0  4 


AMium  cf  «»  NMknd  Rad  Car  ScMSc  RMMKh  OdskH) 


REPORT  DOCUMENTATION  PAGE 


form  4eprovM 
OMC  Vo  0704-0 IM 


t.  ACINCV  usi  ONLY  ((.t4v«  0/4n«;  j.  MAOKT  OATI 

_ _  May  16,  1994 


4.  TI7LI  ANO  SUtTiTLC 


i.  Rcroar  rvat  amo  oatcs  covcrio 
Technical  R 


S.  FUNOINC  NUMBERS 


Electron  Propagator  Theory  and  Application 


i.  AurHOR(S) 


R.  Longo,  B.  Champagne,  and  Y.  Ohrn 


NOOO 14-93- 1-0122 
RT4131072 


7.  FERFORMINO  ORfiANIZATION  NAME<S)  ANO  AOORESS(E$) 

University  of  Florida 
Quantum  Theory  Project,  WM  362 
Gainesville,  FL  32611-8435 


B.  FERFORMING  ORGANIZATION 
RCFORT  NUMBER 


Tech.  Rep.  No.  2 


>.  SFONSORINGi  MONITORING  AGENCY  NAMt(S)  ANO  AOORESS(ES) 

Office  of  Naval  Research 
Chemistry  Program 
800  North  Quincy  St. 

Arlington,  VA  22217-5000 


11.  SUFFICMCNTARY  NOTES 

Submitted  to  Theor.  Chm.  Acta  (1994) 


10.  SFONSORING*  MONITORING 
AGENCY  REFORT  NUMBER 


tZ4.  OiSTRIBUTlON/AVAILABILirY  STATEMENT 

This  document  has  been  approved  for  public  release  and 
sale;  its  distribution  is  unlimited 


12b.  OISTRIBUTION  COOE 


Unlimited 


13.  ABSTRACT  IMwtnum  200  wofttt) 


The  electron  propagator  theory  is  presented  with  somewhat  of  a  historical  perspective  and 
the  working  equations  are  devel^  with  the  aim  to  take  advantage  of  molecular  point  group 
symmetry.  A  new  electron  propagaKR’  code,  tiie  vectorized  electron  propagator  program  OHSP), 
is  introduced  without  full  details  about  its  smictnie  and  capabUities  (such  details  are  being 
published  elsewhere).  Applications  to  the  (UV)  photoelectron  spectra  of  some  donor-acceptor 
complexes  of  boiane  with  carbon  monoxide  and  water  are  presented  at  the  level  of  second-order 
tiieoty  as  an  illustration  of  the  dieoty  and  the  VEP  code. 


17.  SECURITY  CUSSIFICATION 
OF  REFORT 
Unclassified 

NSN  7S40-0'-Za0-SS00 


IB.  SECURITY  CUSSIFICATION 
OF  THIS  FACE 

Unclassified 


It.  SECURITY  CUSSIFICATiON 

o*  abstract 
Unclassified 


IS.  NUMBER  OF  FAGES 


IB.  FRiCE  COOE 


1 20.  limitation  of  ABSTR^ 


StAndara  form  298  R**  2-8 
•nbbchoM  W9  AHV  ^10  fit*  I 

M.IM 


1.  Introduction 


Pn^gators  gained  early  prominence  in  formal  many-body  theory  oi  fermum  systems  (see 
e.g.  [1]  and  references  tlttiein).  Concerns  about  the  eliminati<m  of  unlinked  terms  in  perturbation 
e]q)ansions  and  the  associated  conect  scaling  with  system  size  naturally  led  to  the  propagator 
concept  'neatment  of  double-time  Green’s  functions  [2]  or,  equivalently,  propagators  established 
that  they  provide  a  useful  link  between  quantum  mechanical  treatments  of  pure  state  systems  at 
the  absolute  zero  aiul  that  of  ensembles  at  finite  temperature.  Condensed  matter  theory  employed 
the  propagator  crmcept  to  great  advantage  (see  e.g.  the  review  by  Hedin  and  Lundqvist  [3])  and 
propagator  theory  for  finite  systems  [4,  5,  6]  led  to  new  ideas  for  the  treatment  of  molecular 
systems.  A  detailed  treatment  by  Linderbeig  and  Ohm  [7]  of  a  variety  of  spectrosct^es  for 
molecular  systems  further  developed  and  applied  propagators.  This  paper  focuses  oa  the  electron 
pr(q>agator  theory  and  its  application  to  molecular  systems. 

The  electron  propagator  is  naturally  described  in  Fock  space  in  terms  of  field  qperators. 
Dq)ending  on  only  two  electron  position  and  spin  coordinates  and  a  time  or  energy  parameter 
it  is  in  many  ways  the  t^timal  theoretical  quantity  with  a  dynamical  equatimi  fimn  which  it 
can  be  determined  in  various  approximations.  It  ctmtains  a  wealth  of  informatirm.  Not  rmly 
vertical  electron  binding  energws  (ionizatimi  potentials  and  electron  affinities)  [8,  9]  [10],  but 
also  electron  scattering  amplitudes  [11],  photoionization  intensities  [9,  12],  total  enogies  [13, 
14],  <Hie<]ectron  reduced  density  matrices  [7,  IS]  and  therefore  permitting  the  calculation  of 
cme-dectron  properties  [16].  Derivatives  of  molecular  electrcm  binding  energies  witii  respect  to 
nuclear  displacements  [17, 18]  have  also  been  determined  within  the  electron  propagator  therny. 
Basically  idl  properties  of  a  molecule  are  contained  in  the  propagator.  Thus  electronic  and 
rovibrational  spectra  [10]  and  the  particulars  of  chemical  bon^g  can  be  treated  in  terms  of  tiie 
dectrcm  propagator  [19,  20]. 

Given  a  set  of  orth(»ormal  spin  orbitals  {^(0}  ^  associated  set  of  electron  field 

operators  {ap,<xp}  satisfying  the  anticommutation  relations  (at  equal  times) 

(ap,o,l+  =  [aJ,a{j+  =  [op,a}l+  -  =  0 

tiie  electron  propagator  matrix  is  defined  witii  elements 

«<i,(0;«,(O))  “  -»(« -  ‘’XOM'WWIO) 

Here 

t 

m  =  J  «(T)dT 
—00 

is  the  Heaviside  step  function  expressed  in  terms  of  the  Dirac  delta  fimction, 

N-dectrcm  ground  state.  Atomic  units  are  used  tiuoughout  so  for  instance  h 
Ifeiienberg  equation  of  motkm  for  the  field  q)enton  is 
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where  the  commutator  on  the  ri^t  contains  the  many-electron  HamilttHiian  e^qnessed  in  the 
basis  electron  fidd  (q)eraton 

^M®!***  + 1  S  (M||rj)aJa}a,Or.  (5) 

Pit  ^  .• 


The  one-etectrcm  integrals  involving  the  election  kinetic  energy  and  the  electnm-nuclear  attraction 
terms  are 


and  the  antisymmetric  two-electrcm  integrals  are  {p9||rj)  s  (p9|ra)  —  (p9|ar)  with 

<K|m)  =  /dfl  / <«2#;(£l)^;«2)i^^#r(fl)*.(£j). 

The  electrcm  prc^gator  is  most  commonly  studied  via  its  Fbuiier  transform 

00 

««»;  flj)>r  =  J  «ap(0;  «*(^')))  -  <')» 


(6) 


(7) 


(8) 


which  has  the  spectral  representation  [21,  7] 

“  ,l!5o  ^  If;  +  E„iN)  -  E„,iN  + 1)  +  if, 


<0|«i|m){mloy|0) 


E-Eo{N)^E„,{N 


10) _ 1 

-I)-..,]- 


(9) 


Showing  the  typical  pole  structure,  when  dte  energres  are  discrete  as  they  always  are  in  a  finite 
basis  set  The  energy  dgenstates  are  used  in  this  formal  expression,  i.e. 


fflO) «  £b(A^)|0), 
/r|m>  =  i?„.(JVdbl)|m) 


(10) 


and  only  die  AT-fl  and  A^-1  electron  states  are  involved  as  intermediates  in  die  spectral  repre- 
sratatum. 

It  is  obvious  from  the  spectral  representation  diat  the  electron  {nopagator  has  a  special 
significance  fcr  photoelectron  spectroscopy  and  other  processes  whoe  electron  binding  energies 
are  measured.  The  numerator  contains  the  so  called  Feynman-Dystm  amplitudes 

/,(m)  =  (0|o,|m), 

9pim)  -  <w»|op|0>, 

which  are  obtained  as  residues  at  the  particular  poles  oi  intoest  These  amplitudes  are  important 
in  the  dieoredcal  determination  of  transition  probabilities  for  electrcm  attachment  or  detachment 
processes.  For  instance,  the  intensity  of  a  particular  structure  corresponding  to  the  final  positive 
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ioa  state  m  in  a  photoefectrem  spectrum  of  a  neutral  species  is  proportkmal  to  the  pede  strength 
[8.  9] 

f 

The  propagator  in  equation  (9)  satisfies  fiie  equation  of  motion 

=  (OlK-aJUIO)  +  {{lc„H].-,a\))E 

=  (oilof.  ojl+io)  +  (K;  1^.  “S-Xe 

which  can  be  readily  shown  by  using  die  klentity 

£(£-*r*  =  H-x(i;-xr^  (13) 

with  *  as  Efn(N  + 1)  -  Eo{N)  orx=si  Eq(N)  -  Em(J^  -  1)  in  the  spectral  rqnesentatkm  and 
utiliaing  dm  prop^es  (10)  of  the  energy  eigenstates.  The  two  equivalent  forms  of  the  equadtm 
of  motion  (12)  are  useful  in  manipulating  dm  chain  of  equadems  that  obviously  results  when  dm 
corresponding  equation  of  motion  is  written  for  the  propagator  on  dm  right,  which  involves  dm 
ctmimutator  with  the  Hamiltcmian,  yielding  an  even  more  involved  prq^agatOT  cm  dm  right  with, 
say,  a  double  ctxnmutator  {{a,,  iff]-,  If]-,  and  so  on.  Approximate  treatments  of  dm  propagator 
equations  include  dm  termiiutitm  of  this  chain  of  equatitms  at  snne  suitable  level  or  assumptions 
that  permit  the  summatitm  of  dm  eiqiansions.  This  so  called  decoiqiling  problem  was  addressed 
in  its  generality  by  Un^rberg  and  Ohm  [Q  by  linearuing  the  equation  of  motion.  There  has 
also  qjpeared  discustions  of  truncatitm  schemes  expressed  in  terms  of  diagrammatic  expansitms 
by  Cederbaum  and  Dmncke  [10],  the  moment  conserving  decoupling  using  Padd  qiproximants 
as  well  as  various  renormalization  schenms  discussed  by  Ohm  and  Bom  [22],  and  decoupling 
procedures  expressed  in  terms  of  superoperaiors  [23]. 

Although  equivalent  to  other  procedures  dm  superoperator  formulatkm  provides  a  particularly 
attractive  shorthand  notatitxi.  It  proceeds  by  introducing  a  linear  space  L  of  fennitxi>like  field 
opemtxs 

L  =  {«}, flJ«}or(p  <  ?), oJoJo}o,a,(p  <  g  <  r; <  <  a),  •  •  •},  (14) 

which  supports  a  scalar  product 

(jf|K)  =  (0|ix',y)+|0).  X,Y  e  l.  (is) 

The  superoperator  identity  7  and  superoperator  Hamiltonian  H  are  defined  on  L  such  that 

AX  =  lH,Xl..  ^  ^ 

Iterating  dm  second  form  of  dm  equation  (12)  yields 

-  E  ^  P  ^  ^  (17) 

-(a{|(£/-A)->a}), 

U.  a  matrix  demoit  of  dm  superoperator  resolvent  The  full  matrix  can  be  mqnessed  as 

G{E) «  (a^|(£/  -  (18) 

widi  dm  field  operaton  arranged  in  a  suitable  row  array  on  dm  right  and  cedumn  array  on  dm  left 
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Using  an  inner  jffogection  [24]  manifold  h  €  X  the  superopennor  inverse  can  be  transfonned 
to  a  matrix  invene 

G(i:)  =  (.*|h)(h|(f:/  -  A)hr'(h|a«).  (19) 

This  mqinession  is  the  starting  point  for  the  i^iinaximate  treatmmits  oH  the  electrcm  propagator. 


2.  Approximations  to  the  Electron  Propagator 


Partitioning  of  inner  projection  manifold 

The  simplest  decoupling  of  equation  (12)  is  to  consider  the  so  called  moment  expansion  in 
equadcm  (17)  and  make  the  assumpdcm  that  higher  moments  are  powera  of  die  first  moment 

P  =  (20) 

This  is  (tflimi  lefened  to  as  the  geometric  appraximaticm.  Denoting  (a^|a^)  =  S  the  moment 
expansitm  (17)  becomes 

G(£;) «  £-*sii  +  i;-»s-*p + £;-*s-*FS-*p  +  •  •  •] 

*  S(£;s  -  F)-»S  =  (a*|(JE:7  -  .»)a^)-'  =  Go{E)  ^  ^ 

The  anticommutation  relatkxis  (1)  yielding 

=  =  =  (22) 


(«}!*«♦)  -  <0|(<.„(^,«{J-I+|0) 

have  been  employed,  where  in  die  last  expressimi  is  a  matrix  element  of  die  Fock  operator 
in  the  qiin  oibital  basis.  The  single-particle  reduced  density  matrix  has  die  elmnmits 

=  (01o;o.l0)  (24) 


and  can  be  diagonalized  simultaneously  widi  F  to  give  occupation  numbers  (nr),  U.  7ar  - 
(nr)^$r  and  x^Fx  «=  c  for  some  suitable  unitary  transformation  x  Comparing  equadtxi  (21)  with 
die  formal  qiectral  tqnesentatkm  (9)  it  is  pcnsible  to  write 


Go„{E) 


<»r) 

E  —  tr  — 


»-K)  ij. 


(25) 


For  the  case  diat  die  occiqiation  numben  are  0  or  1  die  rtfmmice  state  |0)  must  be  a  single 
determinant  in  terms  of  die  self-consistent  field  (SCF)  qnn  orbitals 

(26) 

p 
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or  equivalently  eiqiressed  in  terms  of  the  field  operaton  and  die  vacuum  state  as 


|0)  *  JJoJluoc)  s  Ili^lvac).  (27) 

&i  die  following  die  electron  field  operators  {ot  s  r}  refer  to  the  SCF  qnn  oshitals  {xr}.  It 
then  follows  through  (25)  that  at  this  level  of  approximation,  the  occupied  qiin  ortiital  energies 
must  be  identified  with  ionizatkm  pmentials  and  the  unoccupied  ones  with  electrcm  affinities  as 
is  done  diiough  Koc^mans  theorem. 

This  is  omsidered  the  lowest  level  of  approximatirm  and  used  as  a  starting  point  for  all 
hig^  level  treatments  of  the  electron  prc^iagator.  The  S(IF  spin  orbitals  and  then  also  the 
corresponding  electron  field  operators  ns^urally  sqparate  into  an  occupied  set  labeled  by  ... 
and  an  unocciqiied  set  labeled  by  ajb,c, ....  The  labels  p,q.r, ...  refer  to  either  set 

The  inner  projection  manifold  h  used  in  die  expression  (19)  need  (mly  ccmtain  fermimilike 
operaton  (25]  i.e. 

{h}  =  {hi}  U  {ha}  U  {hs}  U  •  •  • 

*  {a^,i^}  U  {a^6^»,»V^o}  U  {o^i^c^tj,t^j^]b^a5}  U  •  •  • 

and  truncadmis  of  this  manifold  corresponds  to  various  iqiproximatimi  schemes.  It  is  ctmvenient 
to  use  an  orthonormal  set  of  inner  projection  basis  elements  [25]  so  that  (hi|hi)  s  i,  and 
(hilh,)  s  0  for  t  ^  i. 

A  fint  step  in  seeking  adequate  approxiroation  schemes  for  die  electnm  propagamr  is 
a  partitioning  of  die  inner  projection  manifold.  When  die  aim  is  to  obtain  a  theoretical 
photoelectnm  spectrum  it  is  oonvmuent  to  choose  the  partitioning 


h*{hi,f} 
f  «  {ha}  U  {ha}  U 


sudi  diat  equation  (19)  becomes 

The  partitioned  form  of  the  inverse  matrix  yields 

G-‘(B)  =  Et-  («»|*a*)  -  (.*|*f)l£l  -  (f  |*f)r‘(fl*«’) 

=  GJ’(£:) -!(£), 

where  die  unperturbed  propagator  and  the  self-energy  tenn  have  been  d^ned  to  show  the  relation 
to  die  so  called  Dyson-like  equatxm  for  die  inopagator[l].  An  untruncated  manifold  f  means  no 
qqnoximatitm  only  a  reformulation  ot  die  [Ropagator  equations.  In  ordm*  to  arrive  at  a  definite 
qiproxiination  and  provide  algoritiuns  for  the  calculatirm  rrf  the  matrix  elements  defining  die 
propagator,  a  reference  state  and  a  truncation  of  die  inner  projection  field  operator  manifold 
must  be  chosen. 


M«cb22,199« 
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Choice  of  reference  state 

The  Hartree-Fock  or  SCP  single  detenninantal  state  |0)  =  \HF)  provides  the  natural  starting 
point  for  die  choioe  of  reference  state  at  varknis  levels  of  i^iproximation.  Althoo^  die  electron 
propagator  dieoty  can  he  fully  develqied  within  a  perturbation  theoiy  framework  with  tmly 
die  SCF  single  detenninantal  reference  state  and  choice  of  inner  projecdtm  manifold  [22]  it  is 
dearable  to  keep  a  balance  between  the  level  of  description  of  the  leferenoe  state  |0)  and  that  of 
die  inner  projecdon  manifold  f.  Such  a  balance  kept  through  various  orders  of  perturbation  theory 
guarantees  hermidcity  the  superoperator  Hamiltonian  matrix  and  the  elimination  of  spurious 
terms  [25].  Starting  frtmi  a  partitioning  of  the  Hamiltonian  and  dius  also  of  the  superoperator 


Hamiltonian 


H  ^  Ho +  SH 


^0  =  y'«i»pV 

»  (32) 


the  reference  state  can  be  expressed  in  terms  of  Rayleigh-Schnddinger  penuibation  theoiy  (RSFT) 
[26,  27],  or  coupled-cluster  (CQ  theoiy.  Also  multiconfigurational  SCF  (MCSCF)  theoiy  has 
been  implemented  [28]  for  the  electron  propagatra*  reference  state.  This  treatise  employs  refuence 
states  based  on  RSPT  and  CC  theory. 

The  hermitici^  problem  consists  of  the  equation 


(X\AY)^iY\AX)* 


(33) 


not  being  satisfied  for  an  iQqnaximate  refeirace  state  or  more  generally  for  a  density  operator 

p  *  po  +  +  •  •  •  +  ^p  (34) 


correa  throngb  order  n  in  perturbation  tberay.  The  average  defining  the  propagator  matrices 
is  dien  a  trace,  ije. 


• 

• 

• 

III 

• 

• 

» 

« 

(35) 

Since 

and  rince 

(36) 

[po,ifo]-  «  0 

(37) 

and 

are  assumed  to  hidd  for  ib  »  1,2, . . .  ,n,  the  error  term  is  of  order  ih-I,  U. 

(38) 

wiry)  -  (y\Ax)’  =  ■n{|{*/,,«H).(x»,yi+}. 

(39) 
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The  Rayleish-SchrOdinger  perturbation  expansim  for  the  reference  state, 

|0)ii5Pr  —  (1  +  /fi  +  ifa  +  ifs  +  •  •  (40) 


is  defined  by 


and 


with  ,  in  particular 


and 


t  • 


»>y  «>s 

=  E  E  . 

i>j>k  m>k>e 


dT’ 


ib?  =  - 
•  2 


lSc  _  \»oiu<g/ 

m  Z-#  m  i* 


(iilljfc); 


iw 


The  denraninators  are  defined  as 


(41) 

(42) 

(43) 


(44) 

(45) 


Df  s=  «i  -  e,  , 


(46) 


and 


i>jf  =  ei  +  e>  -  e,  -  e>  . 


(47) 


The  concept  of  order  in  the  perturbaticm  expansicm  cS  die  electron  prcqiagator  ultimately  means 
order  in  terms  of  the  electmi-electron  interaction  or  equivalendy  two-electron  integrals.  The 
inclusioQ  of  electron  correlation  through  fint  order  in  the  reference  state  is  achieved  with  the 
double  excitation  terms  Kj  whereas  also  the  Ki  terms  are  needed  for  sectmd-order  corrections. 


Coupled-cluster  ratormalization  of  the  retorace  state 

The  coupled-cluster  (CQ  expression  [29]  for  die  reference  state 

|0)cc  -=  ,  (48) 


is  defined  by 


r  =  ri+r2-i-r3-i-...-frjr , 

widi  N  being  the  number  of  elections  the  system,  and 

i  • 


(49) 

(50) 
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(31) 


i>j  •>h 

T,=  Y!,  .  (52) 

i>i>k  a>k>e 

and  90  00. 

Tlie  T  amplitudes  contain  infinite  aider  contiibutitms  to  each  excitation  level  frcmi  tiie  HF 
state.  Thus,  iqplacing  tiie  RSFT IC  amplitudes  with  the  CC  7  amplitudes  may  be  ctmsidered  a 
lenoimaliation  procedure  since  certain  classes  of  penuibation  terms  or  diagrams  are  summed 
to  infinite  onier  [30.  31].  This  idea  has  already  been  employed  in  work  cm  the  polaiizatimi 
propagator  [32,  33]. 

Dctafla  of  the  ooupled-daster  approach 

Expansion  of  die  exponential  in  die  CC  mediod  gives, 

|0)cc  -  (1  +  r,  +  (T,  +  jT?)  + . . 

i  •  »>i  ■>» 

+ 5(E  E  *f“*')(E  E  ) + •  ■  ■ 

»  •  it 

where  die  last  term  can  be  rewritten  as 

5(EE‘f«*‘)(EE‘i**2) = EKE*f<-'o(E‘{»*2)i 

»  •  it  *>i  •  t 

= E  Eiw*)«’*‘’2) + ■ 

»>i  •>t 

and  since  for  the  orthononnal  basis  of  SCF  sinn  orbitals, 

bUa^j  SB  —aUh^j  ,  (55) 

die  CC  reference  state  becomes, 

|o)cc  -  [1 + 7-1  +  (r, +57?)+.. 

=  li+EE‘f<‘*<+EE’-'««*’»*7+"llJ^J^' 

t  ■  »>i  m>k 

where, 

+  «»<}- ^<5.  (57) 

Restricting  die  coi^led-claster  eiqiansion  to  single  and  douUe  ((XSD)  mccitations,  U. 

|0)co  «  WccSD  -  ,  (58) 

die  equations  for  die  T  amplitudes  become  [34], 
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Ti  equation: 


where, 


and. 


-  E‘*0“ii“>  -  5E‘«W“iite)  -  sE'SwiM  • 


Ac 

Tj  equation: 


= wiM) + PM  E  c  (pfc  -  5  E  + 

-  PW)E‘S(^*i  +  5E‘5^*‘)  +  5E’^Jf'^‘“i+ 

+ 5  E  + PW)P(“‘)E  («’*'“«>  -  <:'SMi«»+ 

^  erf  *c 

+ pw)  E‘:(“‘ii'=>)  -  '’ME‘S(“ii‘j>  ■ 

e  Jb 

p- = E‘‘*<*“ii‘*)  -  5  E^“<*'ii‘‘>  ■ 

fee  */c 

p/i  “  E‘S«*ii“> + sE^"  . 

•k  M 

Pj» = . 

M 

Wiii,  =  (WIN;) +P(y)E‘*(‘'ll«>  +  7E«*'II“‘)  ■ 

•  •» 

Wm  -  («‘ll«0  -  P(<|»)  E  (“IM  +  j  E »(?(« IM , 


wmc,  -  (wiicj> +E‘‘(*‘ii“>  -E‘5<*‘i>'»+ 

■  i 


The  ^Eective  two-particle  excitation  opoiattxrs  r  and  f  are 

and, 

—uj  +  2\*  '  ***^/  ’ 


(59) 


(60) 


(61) 

(62) 

(63) 

(64) 

(65) 


(66) 


(67) 

(68) 
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respectively.  In  dw  above  equations,  P(ij)  and  P(ab)  denote  the  permutadonal  interchange 
qperations. 

The  k  coefficients  in  the  RSPT  expulsion  can  be  obtained  from  the  couple-cluster  ansatz 
thnnigh  iteratitn  of  die  T  equations.  For  instance,  the  first  iteration  of  the  T2  equation  gives 
K2,  that  is. 


(69) 


which,  when  applied  in  the  7i  equation  yields 


if=<f(l)  =  5 


Like  i  jkh 


ao) 


Con.«equendy,  it  is  possible  to  write  reference  states  for  the  electrem  propagator  approach  as 
expansion  coefficients  of  the  perturbatitm  therny  or  as  convuged  T  amplitudes  from  die  solution 
of  the  couple-cluster  equations. 

Also,  in  comparing  the  RSPT  and  CC  wave  functitms,  it  is  clear  that 


\0)jiSPT  =  {l  +  Ki  +  K2  +  ...)\HF) 

$ 

\0)ccsD  =  (1  +  It  +  {T2  +  T}/2)  + . .  .)1HF) 

CC-EP  RSPT-^EP 

it 

it  ^  kf 


(71) 


From  the  definititm  of  the  spectral  representation  (9)  it  fdlows  that  the  elements  of  the 
electron  propagator  matrix  G(£)  become  infinite  wten  E  equals  an  electron  binding  energy. 
Then,  the  dements  ci  die  inverse  G'H£)  vanish  at  such  an  energy.  This  result  can  be  used  to 
device  iterative  methods  to  find  the  electron  propagator  poles  and  residues  at  a  given  level  of 
poturbation  expansiem. 


Order  analysis 

Ihincation  of  the  inner  projection  operator  manifold  f  and  the  use  of  RSPT  for  the  refinence 
state  based  rai  the  so  call^  Mpller-Plesset  partitioning  of  the  Hamiltonian  as  eiqjressed  in 
equation  (32)  facilitates  an  order  analysis  of  the  electron  propagator.  In  particular,  the  self¬ 
energy  (see  [26])  and  consequendy  properties  as  e.g.  electrem  binding  energies  and  die  one- 
electttm  reduced  density  matrix  can  be  calculated  to  a  given  desired  oidm  of  electrtm  conelation 
or  electron  interaction.  Starting  from  the  inverse  propagator  matrix  as  given  in  equation  (31) 
and  using  a  shordumd  notation  the  following  etqpresrion  is  tibtained 

(72) 
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where 


IU(£)  =  =  E1-  H.. 

R,/ =  (.*!(£/ -A)f)  =  -H./ 

R/.  =  (f|(B/  -  A)a')  =  -H/.  =  -H*, 

=  (f  l{i!/  -  A)f)  =  -  H,/ 

The  RSPT  expression  for  the  reference  state  then  pennits  the  e^quumcm 


r./=-Eh« 


(74) 


where  the  fact  that  »  0  has  been  used.  Together  with  correspcmding  expanaons  for  the  otho- 

matrices  this  gives  an  expression  for  the  invose  electron  propagator  matrix  through  older  n  as 


g-\e)^-r!S  + 


ehs-  [(i:hi?)(ers?)"(i:h:«)]' 

t*l  \»»1  /  \i«0  /  \»*1  / 


The  first-order  propagator  vanishes  and  the  lowest  order  expressimis  are  obtained  by  calculating 
the  various  matrices  to  specific  orders  and  by  choosing  the  operator  manifold  f  as  follows 


n  =  2  ni  s=  l,n2  =  0;f  =  hs, 

=:2,n2~l;f  —  hs,  (76) 

n  =  4  ni  3,n2  ®  2;  f  »  hs  U  hs, 


which  will  be  further  elabmated  in  the  fdlowing  sections. 

Obviously  the  inversion  of  die  very  large  matrix  Kff{E)  is  one  of  the  difficult  problems  that 
has  to  be  addressed.  An  inversion  could  be  pofoimed  by  employing  a  reduced  linear  equation 
(RLE)  scheme  [35]  but  nqndly  becomes  im]nactica]  wiffi  increasing  basis  sets.  A  number  of 
iq;>proximate  treatments  have  been  proposed  [8, 36, 37]  with  varying  success.  The  order  concept 
can  be  preserved  with  the  identity 

n 


R7>(r)s(RW+£R«)-. 


tvl  iml 


whidi  can  be  iterated  and  truncated  [38,  39,  8]. 


Method  of  solution 

The  general  eiqiression  for  the  elmnent 

G„{E)  ^  {a,\{Ei  ^  Hr\) .  08) 

of  the  electron  propagator  matrix  G{E)  is  symmetry  blocked,  where  each  block  is  formed  by 
the  spin  orbital  indices  p  and  q  belcmging  to  the  same  irreducible  representatkxi  of  the  orffitals. 
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In  ■HHiriftn,  the  electron  pfquigator  matrix  G(E)  is  spin  blocked.  Thmefore,  it  is  sufficient  to 
solve  separately  the  electron  propagator  equaticxis  for  each  symmetry  and  spin  block  pq. 

Fdr  a  given  block  aiMl  energy  E  it  is  pt^sible  to  construct  the  matrix, 

W(£)  =  £1  -  G-*(£)  =  £1  -  (R«(£)  - 

=  E1-  ((a*|(£l  -  A)|at)  -  H.fRj}(E)Hi,) 

=  £1  -  {(£1  -  H„  -  H./R7;(£)H*,) 

»H..  +  H./R7j(£)H*„ 

which  allows  the  expiessitm 

(iJE;-w(f;))G(i;)  =  i.  (80) 


This  shows  that  the  diagonalizadon 


U^W(JE;)U  =  A(£)  = 


Xi(E) 


X,(E) 


0  •  •  • 

=  U'(H„+R’,Rji(£)Rj,)U , 


(81) 


is  important,  with  n  being  the  dimension  of  the  symmetry  pq  block,  and  the  eigenvalue 
conesponding  to  the  spin  orbital  of  intmest  (p)  should  be  the  p-A  pole  (£p)  of  Ae  electron 
propagator  matrix.  Hus  eigenvalue  can  be  used  as  the  next  guess  for  an  iterative  search  of  the 
p-A  pole  or  used  to  obtain  a  guess  for  a  Newttm-RaiAstm  procedure.  Since  Ae  derivatives  of 
W(J?)  wiA  respect  to  can  be  evaluated  analytically,  a  Newton-Raphson  procedure  can  be 
efficiently  employed  to  calculate  the  next  guess  for  £,  so  that  usually,  after  3  or  4  iterations  the 
difference  between  tiie  input  E  and  the  ei^nvalue  is  less  titan  10~^  Hartree. 

Layzer  [40]  treated  such  Qn  general  nonhermitian)  dgenvalue  problems.  When 


Ur(JE;)  =  {1^}.  P=l,2,  --.n 


(82) 


is  the  eigenvector  of  W(£)  corresponding  to  the  eigenvalue  Ar(£:)  and  U*(E)  is  the  eigenvector 
of  coneqtooding  to  eigenvalue  A*(£)  tite  expansion 


G„(E) « 


^  U^{E)U*^{E) 
2^  E-Xr{E) 


(83) 


follows.  The  types  of  possAle  soluticms  have  been  Ascussed  by  Csanak  et  al  [11]  and  details 
have  been  explored  by  Purvis  and  Ohm  [8].  The  pole  of  mferest  Er  is  found  when 


Er  =  Ar(Er) 


(84) 
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and  within  a  finite  basis  the  £r  ue  teal  and  discrete.  Elementary  residue  calculus  gives 


=  [{E  -  Er)G„(E)]  =  rrU^(ErW;,(Er), 


yfAnext 


dAr(£)] 

dE  J  E»Er 


is  die  pote  strengdi  introduced  earlier.  The  resulting  e^qnession 

^  xr^TrU^{Er)U„(Er) 

Gfq{E)  -  2^  YZeL 


can  then  be  compared  to  the  spectral  representation  (9)  to  find 

(0|a,|r)  =  TV^Urr{Er{N  + 1)  -  Eo{N)) 
(r|«,|0)  =  -  Er{N  -  1)). 


(85) 

(86) 

(87) 


(88) 


The  Eeynman-Dyson  amplitudes  direcdy  associated  with  the  various  elecmm  binding  energies 
are  then 


/'’^(f)  =  Ex,(OtV(&)rP 


(89) 


in  terms  of  the  cantmical  molecular  SCF  spin  orbitals. 

The  relationship  of  diese  amplitudes  to  the  electrtm  propagator 


- mmo 


9r(0l/r({') 


- 1)  -  ij?l 


(90) 


iE  +  Eo{N)-Er[N-\'l)  +  iri  ’  E Eo{N)  +  Er(N 
defined  in  positttm  and  spin  coodinate  space  becomes  obvious  finom  the  expansion 

(>«.«)  =  Ex»«KW  (91) 


of  the  fundamental  electron  field  operators  ^((,<}  and  their  adjoints  in  die  basis.  The  Fbynman- 
Dyscm  amplitudes  then  are  of  two  lands  assodated  with  electrtm  attachment  processes  and  with 
electnm  detachment  processes,  respectively  (compare  equadm  (11) 

/ttt)  =  =  '£,Xri(Wr^r'^ 

*  *  (92) 

*(«  - 

f  9 
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Soiutten  Miiaiie 


1.  choose  r 

2.  get  block  pq 

3.  get  guess  for  Er 

4.  coostnict  W{E)  «  Hm  +  Hsf^]}iEr)flL 

5.  diagonalize  UW(iV)  ^  A(£r)l 

6.  get  the  eigenvalue  Ar  associated  with  r 

7.  set  eJSw  *  A, 

8.  get  new  guess  O^ewtcm-Raphstm): 


Ei+i^Ei--G'\Ei)/^ 


dG-HE)\ 

dE  ) 


E»Ei 


(93) 


9.  \Ei  -  JSi-il  <  10”*  7 

a.  ?  yes  Pole  *  Ei 

b.  ?  no  go  to  stq)  3 

Order  analysis  of  the  propagator  matrices 

Second-order  electron  propagator 

In  die  following  die  subscript  1  refen  to  die  hi  »  part  of  the  field  operator  manifold 
and  the  subscript  3  to  die  hs  part  and  so  on.  Through  second  order  die  inverse  of  die  electron 
propagator  matrix  dien  becomes 

Gp5(£:)  =.  Rg>  -  Hi?(Bg>(i:)r>(Hg>)',  (94) 


sriiere 

(hSV  w  -  wm 

(HlV)wi.-0««llu) 

=  {E  +  ti-U-  <»)<«<««(, 

The  second-order  self-eneigy  matrix  then  has  die  dements 


•t®*t 


1  (pg|(«j)falM 


adiae  the  factor  of  1/2  comes  horn  the  relaxatkm  of  die  ordered  indices, 


(95) 


(96) 
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Tliird<ordcr  dcdron  propagator 


A  similar  treatment  through  third  order  yields  [38,  41,  26] 

Gp|(E)  =  Gp|(f:) 

+ h8>  -  (9?) 

-  H<‘>(i45*(£)r'(Hg>)»  -  Hg>{i45'(£)r'(H<«)*, 


adiere 


with 


(•J 


'‘•‘=EE*¥*§. 

t  •>» 

and  where  P(ta)  is  die  permutadonal  interchange  operation.  The  second-order  matrices  are 


(99) 


(h{?W  =  +  (»  -  J’(“‘))EW|m«)«fe 

m>ii  m,c 

(Hg>w = + (» -  pm  E(>“ii“)*5s> 

b>e  m,c 


and  die  first-order  diagcmal  terms  are 

Wlltl)  +  faOt|l<«)  -  <aO»l|t«)  -  fjtmVo)  + 


(JOO) 


(101) 


Fourth  order  and  partial  fourth  order  dectron  propagator 


WidxMit  including  the  operator  manifold  hs  die  full  fourdi-fourth  order  propagator  matrix 
can  be  esqiressed  as 

Gf4‘,(£)  =  Gp5(£:)  +  Hi?  -  HlV(Rg>(£))->Hl?* 

-  -  h1?(I^>(£))->h|1>* 

-  h1?(I^>(£))-'h£>(I^>(£0)‘*h1?* 

-Hg>(Rg>(B))-'Ha>(ag>(£:))->Hll»-Hll>(RS>(B))-'Hg>(Rg>(B))->Hl?' 

(102) 

It  is  general^  more  important  to  include  the  contributions  from  die  hs  manifold  bdore  increasing 
die  order  of  die  expansion  [36]  and  one  dierefore  finds  it  justifiable  to  study  the  etectttm 
propagator  dnough  what  has  been  coined  die  **partia]  fourtfa-ordei^  [19,  37],  where  only  the 
terms  formed  from  die  matrices  alrea^  obtained  in  third  order  are  reodned. 
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Density  matrices  and  many-electron  waveftinction 

The  fint-otder  itdooed  drasity  matxix  or  die  one-matrix  can  be  calculated  from  (90)  as 

= (2»ir'  j  mo-, 

c 

where  the  ccmtour  C  consists  of  the  teal  axis  and  a  semicircle  in  the  upper  half  of  the  ctnnplex 


£-plane  (see  [7]).  Residue  calculus  yields 

r 

which  can  be  compared  to  die  equivalmt  definition 

7(00  =  ^  (105) 

in  terms  of  the  ground  state  Metoctnm  wavefuncdcm  Introducing  die  complete  set  of 
{N  —  l)-election  states  4r,  satisfying 

<i“) 

T 

with  the  compound  coordinate  X  «  Then 

tKI^)  =  ♦o«.X)«;(X)»5(«',X')*,(X')4*4X'  (107) 

and  we  can  identify  the  Feynman-Dyson  amplitudes  with  die  so-called  overUq>  amplitudes 

t.(0  =  /*„(£,  X)*;  (X)dX.  (108) 


Obviously,  if  both  wavefimctions  in  (108)  are  single  determinants  differing  in  one  qnn  ortntal, 
that  qnn  oibital  will  correspcMid  to  the  Feynman-Dyson  amplitude,  while  for  correlated  wavefunc- 
titms  die  amplitudes  are  more  general.  In  additxxi  to  being  important  for  describing  intensities 
in  photoelectrcm  spectra  Feynman-Dyscm  amplitudes  are  also  relevant  for  (e,  2e)  eiqpeiiments 
(see  €.g.  [42]). 

Photoionizatioii  intensities 

The  electrcm  propagator  can  be  used  to  calculate  total  energies,  excitatum  energies,  and  one- 
electrmi  jnoperties  in  general.  However,  it  is  perfaiqis  most  useful  in  die  study  of  photoelecttcm 
qieciroscppy.  In  addition  to  die  electron  binding  mieigies  it  can  also  be  used  to  obtain  estimates 
photoelectron  intrasities. 

The  diffemitial  {dioioionizadon  cross  section,  ije.  die  probability  that  the  systems  absorbs 
one  fdiotoo  causing  a  transition  firom  the  ground  state  widi  wavefunc^  |0)  =  to  an  exdted 
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state  |r}  consisting  of  a  bound  ion  state  with  wavefunction  \N  —  l,r}  s  and  an  electnm 
esci^g  into  the  solid  angle  dO  with  wave  vector  if  is 


N 


(109) 


The  vector  potential  of  die  monochromatic  radiation  field  is  A(f)  ss  with  polaiizaticm 

vector  and  frequency  w.  The  appranmate  form  of  final  (antisymmetric)  state  wavefunctitm 

|r)  «  OAsN^f\{is,iit)*T{iuii.  -  iN.x)  s  (110) 


contains  a  idiotoelectroo  amplitude  v.  which  should  be  a  Feynman-Dyson  continuum  amplimde, 
but  for  a  mcdecular  system  one  usually  has  to  settle  for  an  assumed  fonn.  The  antisymmetric 
projector  is 


N~1 

OAS=‘N-'[l-Y,Pkii]  (111) 

with  the  simple  interchange  permutations  Pks  of  electron  labels.  This  form  of  final  wavefunction 
pennits  us  to  write 

('•iEi(?,)'9)io>-  /  »•(?/, 

>-l  •'  (112) 

+  / 

widi 

jv(f) = (w  -  j  «;(2()i(f,)  •  (113) 

The  seccmd  term  vanishes  if  v(j^/,()  is  strragly  orthogonal  to  •  which  can  be  accomplished 
by  making  v(i/,()  orthogonal  to  the  bound  state  basis.  But  even  when  «rong  ordiogmali^ 
does  not  exist  die  second  term  is  small  for  photoelectnm  enea^s  far  from  dneshold. 

Retaining  only  the  first  term  on  die  right  of  equation  (112)  and  averaging  over  all  incident 
photon  directicms  relative  the  fixed  molecular  frame  and  ova-  die  polarization  dhections  assuming 
arandom  ocioitatitm  of  the  molecules  as  in  a  gaseous  sample  and  for  unpolarized  light,  equation 
(109)  becomes 


^  (114) 

Tb  obtain  this  result  also  the  dipole  qiproximation  was  invoked,  »  1  and  mult^liotion 

by  a  frictar  2  done  to  account  for  die  two  possible  qnn  states  of  the  ejected  elecam  The 
choice  of  a  plane  wave,  an  orthopmalized  plane  wave,  or  a  Coulomb  wave  for  v{if,0  has 
been  tried  with  varying  success  dqmiding  on  the  systnn  and  the  photon  aaetgy  [12,  43].  The 


17 


MHCb22.1994 


ofdiQgooalixed  plane  wave  chnce  aeenu  to  wofk  reasonably  well  for  the  detachment  of  an 
dectioo  fiom  a  negative  kn  leaving  a  neutral  A  plane  wave  is  convenient,  but  not 

particularly  good  under  any  conditions.  It  results  in  [22] 

dor  __  kf 
dn  Srcu 

s.f 


1113) 


where  the  Feynman-Dyson  amplitude  has  bera  expressed  as  a  linear  combination  (rf  atomic 
orbitals  (GTO*s)  {d^} 


The  total  cross  sectitm  is  given  by  integrating  over  all  photoelectron  directions 


ffr  = 


(117) 


Sometimes  die  ratios  of  pole  strengths  are  used  to  predict  relative  intensities  of  structures  in 
a  photoelectron  spectrum  and  that  can  work  for  peaks  of  not  too  disparate  phomelectton  energy. 


The  Vectorized  Electron  Propagator  (VEP)  program 

Approximatioos  to  various  onlets  in  pertuibation  theory  of  die  electron  propagator  have 
been  implemmited  in  the  >fectorized  Electnm  Propagator  (VEP)  program.  Poles  and  associated 
pole  sttengdis  are  computed.  This  code  is  designed  to  be  efficient  by  minimizing  the  number 
of  floating  point  opetadons  and  by  exploiting  vector  and  parallel  features  of  modern  hardware. 
This  is  aooomplisM  by  avoiding  tedhindant  calculadons  dtrou^  the  d^nidon  of  appropriate 
intermediates  and  by  using  symmetry  wherever  suitable. 

The  VEP  code  exploits  die  qiin  and  point  group  symmetry  to  block  die  matrices  fipom  which 
die  propagator  matrix  is  built  Direct  (uoduct  deccmposition  (DFD)  is  used  with  a  sdieme  that 
avo^  redundant  ^mtnetry  checks  before  each  contraction.  This  procedure  makes  die  code 
adqited  to  take  advantage  of  parallel  architBCtures. 

Anodier  feature  of  the  VEP  program  is  possibility  of  renormalized  tteatment  of  die  refierence 
state  of  die  propagator  via  die  use  of  OC  amplitudes  as  well  as  die  standard  MBPT  amplitudes. 
This  is  accomplished  by  interfacing  the  VEP  code  with  die  ACESn  program  systent 

The  DFD  scheme 

Efficient  evaluation  of  matrix  products  necessary  for  die  electron  propagator  calculatums 
can  be  exemplified  by  a  contraction 


Q^TxW.  (118) 

The  indioes  of  die  maoices  can  be  divided  into  two  categories  ,  on  die  one  hand  die  **tatget 
itshoes**  *  *>  which  label  the  matrix  and  the  ^common**  indioes  ei,C2«*  *  *  over  whidi  the 
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oontnction  is  perfonned.  The  target  indices  may  further  be  divided  into  diose  which  Q  shares 
with  T  and  diose  udiidi  it  dunes  with  such  dun 

5Z  ^*(r),eH^e,«(w)  (119) 

Each  etenent  ct  7,  and  W  vanishes  unless  the  direct  product  of  the  inedodble  rquesentatioos 
of  all  indices  labeling  die  quantity  in  questkm  contains  the  totally  symmetric  reinesentaticm.  FCr 
abelian  (sub)groups  this  requires  only  the  following  conditions  to  be  fulfilled: 

r(<i(r))  ®  r(<a(r))  0 .  •  •  =  r(<i(iv))  0  vitjiw))  0  • .  • 
r(^i(r))  ®  r(^2(r))  ®  •  •  •  a=  r(ci)  0  r(<^)  0  •  •  •  (120) 

r(c,)  0  r(c2)  0  •  •  • »  m{W))  0  r(<2(H^))  0  •  •  • 

The  details  of  the  computational  strategy  as  to  storage  of  the  data  structures,  the  manipulatioi  in 
core  memoy,  the  functionality  of  various  subroutines,  etc.  are  rqioited  in  a  separate  publkaticm 
[44]. 

3.  Results 

The  fdiotoelection  spectroscopy  (FES)  techniques  to  study  the  electrouc  structure  of 
atoms  and  molecules  including  transient  and  unstable  qiecies.  This  makes  it  a  powofiil  tool  to 
study  reaction  mechanisms  in  gas  phase  and  oi  surfimes.  PES  has  been  used,  for  instance,  to 
elucidate  the  electronic  structure  ci  donor-acceptor  (DA)  complexes. 

As  an  illustratioi  the  electron  propagator  program  VEP  is  used  to  calculate  the  main  peaks 
of  borane  (BH3)  widi  doncts  sudi  as  H2O  and  00.  These  simple  q^lications  are  limited  to  a 
second-cider  treatment  of  die  dectron  pn;^>agator  (EP2). 

Monomeric  borane  has  a  very  short  life  time,  but  is  a  strong  Lewis  acid  and  may  be  stalnlized 
by  fonning  completes  widi  Lewis  bases  [45].  BH3  resembles  a  transition  m^  atom  in  a 
low  oxidation  state,  in  the  sense  that  it  can  fenn  complexes  widi,  say,  carbtm  monoxide  and 
idioqilioroos  trifiuoride  having  negligible  basicity.  It  has  bemi  suggested  diat  the  complexes 
of  BH3  widi  00  or  PIS  are  formed  via  a  ir-type  delocalization  of  the  BH3  e-ortntals  into 
nnooctqned  00  and  PIS  orUtals.  The  stwty  of  ^sterns  such  as  BH3  •  CO  and  BH3  •  H2O 
may  provide  mqierience  as  to  die  reliability  and  suitability  of  die  electron  propagator  as  a  tool 
for  analyring  PES  experiments  on  heterogeneous  catalysis,  for  instance,  GO  chemisorption  and 
reactioos  (metihantd  tyndiesis)  on  low  axidation  state  transition  metal  oxide  surfeces,  such 
u  ZnO(10lO)  and  also  010(111)  [46]. 

Molecular  Geometries 

Hie  electroo  propagamr  program  is  implemented  in  the  ACES  n  program  systnn  [47].  All 
calctdation  presented  were  perfbrmed  on  an  IBM  RS/60(X)-580.  The  baris  sets  are  correlated 
conristent  pVDZ  [48],  which  conrist  of  (9s4pld^3s2pl^  for  first  row  elements  and  (4slp)/[2slp] 
for  hydrogen.  AH  stroctmes  are  optimized  m  die  RHF/MBPT[2]  levd  of  theory.  The  results 
are  Hsied  in  lUde  1. 
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Ikble  1  Optimized  feometriet  at  the  MBPTI2]  ievd  with  pVDZ  bates.  The  nootion 
meant  that  die  pratoa  liet  in  the  tymmetiy  jdane  of  the  Q  point  groiq>. 


Complex 

Pmnt  Group 

Geometrical 

Parameter 

MBPTI21/i)VI>Z 

Experimental 

BC 

1.565  A 

1340  A 

HbB-CO 

C3, 

CO 

BH 

1.144  A 

1.217  A 

1.131  A 

1.194  A 

CBH 

103.8* 

104.2® 

BO 

2385  A 

HiBOC 

C3. 

CO 

BH 

1.147  A 

1.203  A 

OBH 

90.9® 

BO 

1.730  A 

OH 

0.963  A 

HsBOHi 

C. 

BH&BHt 

1.214 

OBH 

100.6® 

OBH( 

103.7® 

R3B 

Dsb 

BH 

1.203  A 

CO 

Coov 

CO 

1.147  A 

1.128  A 

H2O 

OH 

0.964  A 

0.958  A 

HOH 

102.0® 

1043® 

Symmetry 

The  effects  of  motecular  symmetry  oa  Ae  peiAnnance  ttf  the  code  is  Ulus* 
trsted  by  cateulating  Ae  PES  BH3  at  the  EP2  level  in  Ae  pVIZ  basis,  Le. 
(10sSp2dlfV[4s3p2dlf|//(Ss2pldV[3s2pld]  m  Ae  maximal  Abelian  subgroup  C2v  and  m 
Cl.  The  Aeoietical  factors  (rf  ledncdon  (in  computing  time  or  total  number  of  floating  point 
operations)  doe  to  symmetry  (FRS)  have  been  defined  [34]  and  is  found  in  many  cases  for  total 
energy  calculations  using  flto  OCSD  (coo|ded-cluster  singles  and  doubles)  levd  of  Aeory  to  be 
dose  to  die  symmetry  gtoiqt  order  square.  As  diqilayed  in  Thble  2  Ae  presoit  calculation  does 
not  quite  leadi  fliat  efficioicy  impiovemait  wiA  Ae  symmetry  treatment,  but  sdU  a  respectable 
fedncdon  in  computing  tiine  is  achieved. 
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lUilB  2  An  FRS  It  a  "Uctpi  of  reducdoo  (be  to  symmetry”  of  cpu  time  used.  The  VEP  code  b  the  ACES  n 
pngnm  lysiein  adiieves  fot  tbs  case  a  facttr  of  more  than  10  qwed  q>  b  the  executbo  oi  the  code. 


Molecule 

Sire  basis 

Pennt  (Sroop 

Achieved  FRS 

BH3 

72 

Cl 

1.0 

C2v 

10.3 

Donor-Acceptor  PES 

The  photoelectrra  qiectnim  with  a  UV  source  (UPS)  of  H3B*OH2  calculated  at  the  SCF 
(Koopmans*  theorem)  and  tiie  EP2  levels  are  cmnpared  to  experimental  results  [49]  in  Thble  3. 

Ihble  3  Electron  propagator  at  the  SCF  level  (Koopmana)  and  at 

the  EP2  level  are  conquied  widi  cxperbieat  for  the  H3BOH2  comidex. 

Assignment 

Koopmans  (eV) 

EP2(cV) 

(Pde  Strength) 

Experiement  UPS 
(eV) 

ir(B-H)  7a' 

11.3 

10.5  (0.93) 

9.7 

x(B-H)7a" 

11.8 

11.1  (0.94) 

10.6 

<r(B43)6a' 

15.6 

133  (0.91) 

11.8 

n(0)  Sa' 

16.4 

14.2  (0.90) 

13.2 

n(0)  +  <7(B0)  4a' 

20.7 

18.6  (0.89) 

14.4 

Due  to  the  hydroly^  of  diborane  the  expaimental  investigation  of  tiie  H3BOH2  complex 
is  difficult  [45]  and  introduces  some  uncertainty  about  whether  tiie  observed  features  in  the  Hel 
qtectrum  really  is  due  to  lisBOHi  or  smnetfaing  else.  The  agreement  between  tiie  calculated 
(EP2)  peaks  and  tiie  UPS  qiectrum  is  as  expected  excqit  for  the  observed  feature  at  14.4  eV. 
This  is  not  con^tent  with  the  theoretical  r^t,  but  before  suggesting  tiiat  this  feature  might 
not  be  due  to  H3BOH2  the  dectron  propagator  calculations  have  to  be  carried  to  the  third  or 
partial  fourth  order  and  also  a  larger  basis  used. 

Comparison  of  tiie  Koopmans*  theorem  and  the  EP2  results  with  experiment  for  tiie  HsB^CO 
complex  is  presented  in  Tl^  4. 
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Ikble  4  Ekcttoo  prapigMar  poka  m.  tfae  SC7  (Kaqmns)  levd  and  at 
diB  EP2  levd  ooiqMnd  with  otperimeDt  for  tbe  1%B<X)  eomidcat. 


Assignment 

Koopmans’  (eV) 

EP2(eV) 
(Pole  Strentfa) 

Experiment  UPS  (eV) 

2e(B-H) 

12.8 

11.9  (0.92) 

11.9 

6ai  (B-Q 

15.1 

13.9  (0.91) 

14.1 

le  (C^) 

18.6 

17.1  (0.87) 

17.0 

5ai  (GO) 

21.6 

18.1  (0.85) 

18J* 

«  Adiabatic  ionizatioa  enosy 


Similar  calculations  are  carried  out  for  the  HaBOC  complex  and  are  reported  in  Ikble  5. 


Ikble  S  Compariioo  of  the  dectron  propagMor  poles  at  tfae  SCP  goopeaaM) 
levd  md  at  the  EP2  levd  with  the  same  UPS  qMCtnan  as  in  IkUe  4. 


Assignment 

Koopmans  (eV) 

EP2(eV) 

(Pole  strength) 

Experiment  UPS  (eV) 

2e(B-H) 

13.2 

12.6  (a94) 

11.9 

6ai  (B-Q 

’  153 

14.0  (a92) 

14.1 

Ic  (C-O) 

173 

163  (0.89) 

17.0 

5ai  (B-H) 

18.7 

173  (0.92) 

183* 

4ai  (C-O) 

22J 

19.0  (0.86) 

^Adiabatic  ioaizatkm  energy 


There  is  definitely  a  better  agreenrent  betwem  tfae  calculated  and  the  observed  PES  for 
HsB^GO  dian  for  HsB^OC  indicating  die  discriminatory  power  of  die  electnm  propagator  dieory 
even  at  this  primitive  level 

In  coBclurion  one  can  again  reafifirm  v^utt  already  has  been  established  by  many  workers 
in  die  field,  namely  diat  die  propagator  dieory  is  an  ^ipropxiaie  and  practical  i^iproach  to  the 
interpretatun  and  prediction  qiectn.  Ihe  tesulo  presented  here  also  show  that  in  odor 
to  contain  truly  quantitative  agreemoit  widi  eiqperimmit  it  is  necessary  to  consider  electron 
propagator  dieory  at  die  dtird  and  partial  fourdi  order  and  to  also  be  able  to  accommodate  larger 
batis  sets. 
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Tables 


lUde  1  Opdmized  seomraiei  at  tbe  MBPT[2]  levd  with  pVDZ  tmu.  The  notttion 
B«  ineani  that  the  pratoo  lies  in  the  lymmetiy  ptaae  of  the  €•  point  grenp. 


Qmiplex 

Pdnt  Gronp 

Geometrical 

Parameter 

MBFr[2]/pVDZ 

Ejqterimmital 

BC 

U65  A 

1J40A 

HsB^CO 

C3v 

CO 

BH 

1.144  A 

1.217  A 

1.131  A 

1.194  A 

CBH 

103.8* 

104J1* 

BO 

2.585  A 

H3BOC 

C3r 

CO 

BH 

1.147  A 

1.203  A 

OBH 

90.9“ 

BO 

1.730  A 

OH 

0.963  A 

H3BOH2 

c, 

BH&BH, 

1.214 

OBH 

100.6* 

OBH, 

103.r 

H3B 

Dst 

BH 

1J03  A 

CO 

Coov 

CO 

1.147  A 

1.128  A 

H2O 

C2T 

OH 

0.964  A 

0.958  A 

HOH 

102.0“ 

104.5* 

lkUe2  AnFRSiia‘^Ktarafitdiiction<laeiDiyiiUBi»iy”ofcpatiiiieaMd.  Hie  VEP  code  in  die  ACES  n 
procmn  system  achieves  fbr  dus  case  a  bctor  of  more  than  10  qieed  op  in  the  execution  of  the  code. 


Mdecule 

Size  of  basis 

Point  Group 

Achieved  FRS 

BHa 

72 

Cl 

1.0 

103 
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Ikble  3  Ekctioo  prapagMor  pokt  tt  die  SCF  level  (Koopmeni)  nd  at 
tbe  EP2  levd  are  coiniiared  with  experiment  liar  the  HsBKMIi  complex. 


Assignment 

Kot^amans  (eV) 

EP2(eV) 

(Pole  Strength) 

Eacpetiement  UPS 
(cV) 

ir(B-H)  7tl 

11.3 

lOJ  (0.93) 

9.7 

7a" 

11.8 

11.1  (a94) 

10.6 

<T(B-0)6a' 

1S.6 

13J  (0.91) 

11.8 

n(0)  5ti 

16.4 

142  (0.90) 

13.2 

n(0)  +  rr(B-0)  4a' 

20.7 

18.6  (0.89) 

14.4 

4  V  * 


IkUe  4  Ekctran  prapiguor  pota  ic  tlie  SCF  (Kbopnans)  levd  and  at 
ite  EP2  ievtf  compared  with  experimeat  Cor  the  %B*00  compta. 


Assignment 

Koopmans*  (eV) 

EP2(cV) 

(Pol6  Strentfa) 

Experiment  UIK  (eV) 

2e(B>H) 

1Z8 

11.9  (0.92) 

11.9 

6fti  (B-C) 

15.1 

13.9  (0.91) 

14.1 

le  (C-O) 

18.6 

17.1  (0.87) 

17.0 

5ai  (00) 

21.6 

18.1  (0.85) 

18.5* 

*  AdiiJMttic  iauzadon  enetgy 


4 


Mhreh24.19»4 


«  •  ' 


lUrie  S  CompiriMo  of  the  okctran  prapogaior  polei  u  (be  SCF  (KoopmoM) 
level  nd  at  tbe  EP2  level  with  tbe  nme  UPS  ^tecimin  at  in  lUde  4. 


Assignment 

Koqmians  (eV) 

EP2(cV) 

(Pole  strength) 

Ejqpeiiment  UPS  (eV) 

2e(B-H) 

132 

12.6  (0.94) 

11.9 

6ai  (fi-Q 

15.3 

14.0  (a92) 

14.1 

le  (C-0) 

17J 

163  (a89) 

17.0 

Sai  (B-H) 

18.7 

173  (0.92) 

183* 

4ai  (C-0) 

215 

19.0  (0.86) 

♦Adiabatic  ioaizatkn  anagy 


* 


5 


MaRh24.1994 


